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Abstract 

Using three-dimensional particle-in-cell simulations we show that a strongly nonlinear plasma 
wave excited by an ultrahigh intensity laser pulse works as a compact high-brightness source of X- 
ray radiation. It has been recently suggested by A. Pukhov and J. Meyer-ter-Vehn, Appl. Phys. B 
74, 355 (2002), that in a strongly nonlinear regime the plasma wave transforms to a "bubble", 
which is almost free from background electrons. Inside the bubble, a dense bunch of relativistic 
electrons is produced. These accelerated electrons make betatron oscillations in the transverse 
fields of the bubble and emit a bright broadband X-ray radiation with a maximum about 50 keV. 
The emission is confined to a small angle of about 0.1 rad. In addition, we make simulations of 
X-ray generation by an external 28.5-GeV electron bunch injected into the bubble. 7-quanta with 
up to GeV energies are observed in the simulation in a good agreement with analytical results. 
The energy conversion is efficient, leading to a significant stopping of the electron bunch over 5 mm 
interaction distance. 
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The development of novel high-brightness compact X-ray sources is important for many 
research, industrial and medical applications. Synchrotron light sources (SLSs) are the most 
intense X-ray sources today. In an SLS, the radiation is generated as a result of relativistic 

n 

electrons scattering by a bending magnetj_ magnetic undulators or wigglers [1], or by high- 



power laser pulses (Compton scattering) . Recent experiments, which explore the 

interaction of an intense 28.5-GeV electron beam with plasma at Stanford Linear Accelerator 
Ce,.e. (SLAC) 33. have show. *at a„ ion channel can be snecessfnlly used as a wiggle, 
to produce the broadband X-ray radiation: the electron beam propagating in plasma blows 
out the background electrons and generates an ion channel. 

A relativistic electron running along the ion channel undergoes betatron oscillations about 
the channel axis due to the restoring force acting on the electron by the uncompensated ion 
charge. For small amplitudes, rofcb <^ 1, the betatron oscillation is close to a harmonic 
motion with the betatron frequency Ub = ck^ = io^j^/T?). Here tq is the radial excursion 



of the electron, cUp = ^J^n^TL^m is the background plasma frequency, Uf. is the elec- 
tron density that is equal to the ion charge density in the channel, 7 is the relativistic 
Lorentz factor of the electron, e is the electron charge, m is the electron mass and c is 
the speed of light. Relativistic electrons executing betatron oscillations in the ion channel 
emit short-wavelength EM radiation ^, 01 . Some features of this radiation spectrum have 
been studied in the recent publications j?, Q|. The fundamental wavelength of this radi- 
ation is close to A ^ Afe/ (27^) for small-amplitude near-axis betatron oscillations, where 
\h = 2TT/kb. The emission at the fundamental frequency has been call ion channel laser 
in the work of D. Whittum et al If the amplitude of the betatron oscillations be- 

comes large, then the electron radiates high harmonics. If the plasma wiggler strength, 
K = 7fcfero = 1.33 x lO^^^^y^yrif, [cm~'^]ro [/im], is so high that K ^ 1, then the radiation 
spectrum becomes quasi-continuous broadband. It is similar to the synchrotron spectrum, 
which is determined by the universal function S{uj/uJc), where S (x) = x K^/z (0 
ujc is the critical frequency For frequencies well below the critical frequency {uj <^ cUc), 
the spectral intensity increases with frequency as cu^/^, reaches a maximum at ~ 0.29ci;c, and 
then drops exponentially to zero above Uc- The critical frequency for a relativistic electron in 
an ion channel is huJc = (3/2)7^fcroA;^ ~ 5 x 10~^"'^7^ne [cm~^] rg [/im] MeV. The synchrotron 
radiation emitted from an ion channel has been observed in a recent experiment 

The synchrotron radiation is confined to a narrow angle Or ~ K/'j because of the strongly 
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relativistic motion of the electron. The averaged total power radiated by an electron un- 
dergoing betatron oscillations is (Ptotai) — e'^c-f'^kpTQ/ 12. It follows from this expression 
that the radiated power is proportional to the squared density of ions in the channel. This 
fact has been confirmed in the experiment 9]. The averaged number of photons with the 
mean energy HlUc emitted by the electron is {Nx) ^ {2n /9){e^ /nc)NoK ~ 5.6 x IQ-^A^'oi^, 
where A^o is the number of betatron oscillations executed by the electron. 

In the SLAC experiment the ion channel has been produced by the electron beam 
itself in the blow-out regime when the electron beam density, n^, is higher than the 
plasma density. The density of a relativistic electron beam cannot be very high because of 
the technology reasons. This leads to a serious limitation on the gain in the radiated power, 
which is quadratic in plasma density. 

The use of a high-power laser could overcome this limitation. The high-power laser 
pulse can expel plasma electrons by its ponderomotive force and create the ion channel Q| . 
Moreover, the strongly nonlinear broken-wave regime has been recently observed in 3D PIC 
simulations Q]. In this regime, the background electrons are completely evacuated from 
the first period of the plasma wave excited behind the laser pulse, and an "electron bubble" 
is formed. The ion density in this bubble is many orders of magnitude higher than that in a 
simple beam-plasma interaction. For example, the ion density in the laser-produced channel 
can be as high as 10^^ cm~^ jisl . [l4 1. This is 10^ times higher than that in the recently 
reported beam-plasma experiment p. Therefore the radiated power in the laser-produced 
channel may be 10^° times higher. The bubble moves with the group velocity of the laser 
pulse, which is close to the speed of light. A relativistic electron bunch injected into the 
bubble can propagate inside the bubble over a very long distance. Hence, in spite of the 
small length of the bubble itself, the electrons can oscillate in the bubble for a long time. 

It has been recently shown by three-dimensional particle-in-cell (PIC) simulations that 
a dense quasi-monoenergetic bunch of relativistic electrons, collected from the background 

n 

plasma, can be generated inside the bubble [1^. Because of the bubble focusing the bunch 
has a much higher density than the background plasma. In the present work we show that 
betatron oscillations of the bunch in the transverse fields of the bubble lead to the efficient 
X-ray generation, which can be used for the developing of table-top high-brightness X-ray 
radiation sources. 

We perform a numerical simulation of the X-ray generation in laser-plasma interactions for 
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the strongly nonlinear broken-wave regime when the bubble is formed behind the laser pulse. 
We use the fully 3D PIC code Virtual Laser-Plasma Laboratory Q| . The incident laser pulse 
is circularly polarized, has the Gaussian envelope a{t,r) = aoexp(— r^/r| — t'^/Tl), and the 
wavelength A = 0.82 fim. Here a = eA/mc^ is the relativistic laser amplitude, = 8.2 fim, 
Tl = 22 fs, Oq = 10. The laser pulse propagates in a plasma with the density rig = 10^^ cm~^. 

Fig. m presents snapshots of the laser pulse (the colored scale) and the electron density 
(the black/white scale) at different distances. The laser pulse has passed 14 Rayleigh lengths 
{Z^ = 7rr|/A) after the interaction time Tint = 4500A/c. Thus, the lifetime of the bubble is 
about 3500 A/c ~ 10 ps. Electrons, trapped in the bubble, form the relativistic bunch. We 
observe as the bubble stretches and the bunch elongates with time. 

Despite the fact that the bunch density is higher than the background ion density, the 
transverse force acting on the accelerated electrons, F±, is mainly determined by the electro- 
static focusing force from the ions, see Fig.|2K- This is because the charge force of relativistic 
electrons and the self-generated magnetic force almost cancel each other The energy 
spectrum of the electron bunch is shown in Fig. Eb- We observe formation of the quasi- 
monoenergetic peak 3]- At ct = 4000A the peak is located at 360 MeV. We calculate the 
corresponding wiggler strength of ~ 89 ^ 1. Thus, the electrons emit X-rays in the 
synchrotron regime. The number of electrons in the bunch is about 6.5 x 10^° at this time. 
The total energy of electrons of the bunch is about 3.3 J that is about 20% of the laser pulse 
energy. The number of betatron oscillations experienced by the electrons up to this time 
was A^o = cTint/Xb - 8.6. 

To simulate the X-ray generation we suppose that at any given moment of time, the elec- 
tron radiation spectrum is synchrotron-like ^J. The spectrum integrated over solid angle is 
defined by S{uj/lJc)- The critical frequency Uc is given by the relation Uc = (3/2)7^ |Fj_|/(mc), 
F± is the transversal to the electron momentum force. In our PIC code, we follow trajecto- 
ries of each electron and calculate the emission during the interaction. The emitted radiation 
ex.. a .00,. on .he e.ec.on g TKe .eooU fo.ce wa. .eluded ,„to .Ke equation, of e.ec..o„ 
motion in our simulations. 

The synchrotron spectra at ct = 1000 A and ct = 4500 A are presented in Figs. El (a, 
b). The surfaces shown in Figs. El (a,b) give the number of photons within 0.1% of the 
bandwidth [Afiw = lO'^fuj) per solid angle, 27r sin 6'(i6': Nx = Aud"^ Nx /{27t sin 9 dudO). It 
is seen from Fig. O (b) that the relativistic bunch radiates highly energetic photons within a 
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very narrow cone. The maximum of the radiation spectrum is located at about 50 keV. The 
analytical estimates for electron energy predict the maximum of S{x) ~ O.SfkJc — 55 keV 
that is in a good agreement with the numerical simulation data. It is seen from Fig. |21 that 
the radiation from the bunch is confined within the angle 6' ^ 0.1 rad and the theoretical 
estimate is about 0.2 rad. The photon fiux (the number of photons per second in 0.1% 
bandwidth) and the spectral brilliance of the source at ct = lOOOA and ct = 4500A are 
shown in Figs. 01 (c, d). We can estimate the fiux and the brilliance using the following 
formulas 0| $ ~ (Auc/uJc) Nx{c/Lf,) and B ~ $/(47r^6'|S'|) , where Lb is the bunch length, 
S'i? ~ vr [rf + c^Tj^i6'|j/(47r^)] is the effective source size of the radiation and rb is the bunch 
radius. 

To emphasize the advantage of the X-ray generation in the laser-produced ion channel in 
comparison with that in the self-generated channel, we perform a numerical simulation of 
the X-ray emission from an external 28.5-GeV electron bunch. The bunch has a diameter 
2ro = 24.6 fim and a length Lb = 82 fim with the total charge Qb = 5.4 nC. The plasma 
and laser pulse parameters are the same as in the previous simulation. The electron beam 
density was much smaller than that of the background plasma, so that the laser pulse and 
bubble dynamics is not strongly affected by the external electron bunch. At the beginning of 
interaction the front of the electron bunch is close to the center of the laser pulse (see Figs. |3] 
(a)). The head of the bunch has overtaken the laser center by some 46A after the interaction 
time Tint = 4500A/c. The number of betatron oscillations during the interaction time was 
A^o = cTint/Xb — 1.1. It is seen from Fig. HI (b) that the laser pulse and the bubble remain 
structurally stable during the full interaction and the bunch is focused at this moment of 
time. 

The synchrotron spectrum after the interaction time Tint = 4500 A/c is presented in 
Fig. El (a). In the present simulation we do not consider the emission from the background 
plasma electrons. At the given plasma density, the plasma wiggler strength parameter is 
about K ~ 817. It is seen from Fig. El (a) that the relativistic bunch radiates highly energetic 
photons within a very narrow cone. The maximum of the bunch radiation spectrum is located 
at about 210 MeV. The analytical estimates predict the maximum of S{x) ^ O.SHuJc — 385 
MeV. The disagreement is caused by the bunch stopping because of the radiation damping 
force. We also observe a significant photon fiux up to the energy of 10 GeV. The radiation 
from the bunch is confined within the angle ^ ^ 10 mrad that is close to the theoretical 
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estimate 15 mrad. The total number of photons emitted by the bunch are about 2 x 10^^. 
This means that every electron of the bunch emits about 6 photons. The estimation for the 
photon number with the critical frequency is Nx = N^, (Nx), where is the number of 
electrons in the bunch. The estimation is in a good agreement with the numerical simulation 
results. The bunch lost about one third of its energy after Tint- The energy distribution of 
the bunch electrons after the interaction is shown in Fig. El (b). The photon flux and the 
brilliance versus the photon energy are shown in Figs.|Hl(c,d). The brilliance at the beginning 
of interaction is slightly higher than at the end because, at the beginning, the bunch is not 
yet focused and, therefore, emits at small angles. It follows from the Figs. El (c,d) that the 
photon energy, flux and brilliance of the X-ray emission from laser-produced ion channel are 
several orders of magnitude higher than the ones observed in the self-generated ion channel 



In Conclusion, we propose a novel compact and intense x-ray radiation source based on 
the strongly nonlinear broken-wave laser-plasma interaction. The brilliance of this source 
is some two orders of magnitude higher than that of the best x-ray sources available today 
In addition, the radiation is polychromatic, covers that multi-keV range and comes in 
sub- 100 fs pulses. This bright novel source of femtosecond X-ray pulses will have important 
scientific applications by enabling the direct measurement of atomic motion and structural 
dynamics in condensed matter on the fundamental time scale of a vibrational period P|. 
The 100 fs time scale is characteristic for atomic motion associated with ultrafast chemical 
reactions, nonequilibrium phase transitions, surface dynamics and even ultrafast biological 
processes. 

In addition, our proposed radiation source provides a sufficient number of photons per 
pulse to carry out these studies in a single-shot regime. This is crucially important, e.g., for 
the biological processes. The polychromaticity of the radiation source may allow to probe 
simultaneously different atomic species in complex or disordered materials. 

The proposed radiation source can be a table-top laser, and the plasma interaction length 
is less than a centimeter. This must be compared with many 100 m structures of the 



conventional synchrotron sources [l^. The high ion density in the laser-plasma wiggler 
provides several orders of magnitudes higher energies of the x-ray photons than that observed 
in the recent experiment with self-generated ion channels and in the designed FELs. 
One of the authors (I. K.) gratefully acknowledges the hospitality of the Institute for The- 
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FIG. 1: (Colour) The evolution of the laser pulse intensity (the coloured scale) and the bubble (the 
electron density is given in the black/white scale) in the strongly nonlinear broken- wave regime. 
The laser pulse propagates in a plasma layer from left to right. The plasma density and the laser 
intensity at a) ct/\ = 500, b) ct/\ = 1000, c) ct/X = 2000, d) ct/X = 4000. 
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FIG. 2: (Colour) a) The transversal force acting on the relativistic electrons moving in the x- 
direction at ct/X = 2000. b) Temporal variation of the energy spectrum of the electron bunch: (1) 
ct/X = 1000, (2) ct/X = 2000, (3) ct/X = 3000, (4) ct/X = 4000. 
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FIG. 3: a) The synchrotron spectrum from the plasma at ct/\ = 1000, b) at ct/X = 4500, c) 
the photon flux (the number of photons per second in 0.1% bandwidth), d) the spectral brilliance. 
The dashed line in frames c) and d) corresponds to ct/\ = 1000, the solid line corresponds to 
ct/X = 4500. 
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FIG. 4: (Colour) Temporal evolution of the plasma density, laser intensity and the envelope of the 
external 28.5-GeV electron bunch (blue): a) at the beginning of interaction and b) at ct/A = 1500. 
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FIG. 5: a) Synchrotron spectrum from the external 28.5-GeV electron bunch at ct/X = 4500; b) 
energy distribution of the bunch electrons: the solid line corresponds to ct/\ = 4500, the dashed 
arrow marks the initial energy of the electron bunch; c) photon flux and d) spectral brilliance. 
The dashed line in frames c) and d) corresponds to ct/\ = 500, the solid line corresponds to 
ct/X = 4500. 
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